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Abstract 

In this work we revisit the growth of small primordial black holes 
(PBHs) immersed in a quintessential field and/or radiation to the su- 
permassive black hole (SMBHs) scale. We show the difficulties of sce- 
narios in which such huge growth is possible. For that purpose we 
evaluated analytical solutions of the differential equations (describ- 
ing mass evolution) and point out the strong fine tuning for that 
conclusions. The timescale for growth in a model with a constant 
quintessence flux is calculated and we show that it is much bigger 
than the Hubble time. The fractional gain of the mass is further eval- 
uated in other forms, including quintessence and/or radiation. We 
calculate the cosmological density due to quintessence necessary to 
grow BHs to the supermassive range and show it to be much bigger 
than one. We also describe the set of complete equations analyzing 
the evolution of the BH+quintessence universe, showing some inter- 
esting effects such the quenching of the BH mass growth due to the 
evolution of the background energy. Additional constraints obtained 
by using the Holographic Bound are also described. The general equi- 
librium conditions for evaporating/accreting black holes evolving in a 
quintessence/radiation universe are discussed in the Appendix. 
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1 Introduction 



Primordial black holes may be important in a number of situations. 
If PBHs formed in the early universe from primordial perturbations 
or some other mechanism, it is still not excluded that a population 
may survive. Some of them may even become the supermassive black 
holes that energize active galactic nuclei driving the powerful source of 
quasars, and may be a powerful sources of gravitational waves as well 
(see the review of M.Rccs [1]). 

The origin of supermassive black holes (SMBHs) remains a mystery, 
and since there are robust observational evidence for their existence 
with huge masses (10 6 M Q — 1O 8 M ), a big problem arises. How these 
objects can be formed in the universe? Were their initial masses very 
large from the beginning? There arc several proposals for understand- 
ing these questions. One can suppose that the big black holes lurked 
in the AGNs can be formed from a primordial end point of the first 
generation of stars. A galactic nucleus could be populated by several 
massive bodies, and these objects fused into a only and single mas- 
sive collapsed body, forming a central black hole. These models are 
discussed, for example, in Refs.[2] and references therein. 

A second class of models (studied here) assume a huge growth of 
PBHs as a working hypothesis. Quite independently of the formation 
scenarios of small black hole seeds, the important question is whether 
they can grow to the supermassive range by means of accretion. 

The analysis of the Hawking evaporation [3] of PBHs has leaded 
to conclude that PBHs initially larger than M haw <~ 10 15 g could have 
chance to survive and contribute to a small fraction of the dark matter. 
A subset of PBHs may have grown to the supermassive range residing 
within the observed AGNs. In Ref. [4] we evaluated in detail the 
evolution and survival of PBHs taking into account evaporation and 
absorption from a radiation-dominated environment. We proved that 
all those PBHs within in the radiation-dominated era remained colder 
than the environment, and had a very small mass growth. 

The simplest possibility for reaching the supermassive range is that 
PBHs formed in the early universe accreted mass from an energetic 
radiation environment (the Hawking evaporation was important only 
when their proper temperature became higher than the background 
temperature). If PBHs were massive enough, they could be colder 
than the environment, and accretion will drive their evolution. The 
accretion term may be constructed by multiplying the cross section for 
classical absorption of rclativistic particles by the flux coming through 
the horizon. Taking into account the gravitational focusing of the 
incident radiation, this quantity is of order 
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(dM\ 27tt 2 

where the factor 27-k/A describes the relativistic beaming of parti- 
cles, see Zel'dovich and Novikov for details [5]. Since the gravitational 
radius is R g = 2GM/c 2 7 this accretion is proportional to the second 
power of the black hole mass M. The formula above assumes that (in a 
first approximation) classical accretion does not affect the metric that 
describes the black hole itself. The issue is whether one can construct 
quintessence models which not only give the correct dynamics for the 
expansion, but also make PBHs grow to the supermassive range. This 
possibility is the subject of the following Sections. 



2 Quintessence models and Supermassive 
Black Hole growth 

There is now very good evidence that a dark energy is present in our 
universe and represents ~ 0.7 of the critical density. Today, the big 
challenge is understand the physical nature of this contribution. A cru- 
cial ingredient is that, contrary to the dark matter, dark energy must 
possess a negative pressure as required for an accelerating universe. 
Then, g e ff = (g + 3P) < 0, implies to < —1/3, with w = P/g the 
velocity of sound. 

At first sight, the most natural candidate is the cosmological con- 
stant. The cosmological constant can be viewed as a fluid exactly pos- 
sessing to = —1. Since g+3H( 1+ to) g = 0, implies that g = constant. 
However, as it is well-known, explaining the dark energy with a cos- 
mological constant runs into problems. The first, is that we need to 
generate a tiny ("unnatural") value at present epoch, to explain the 
magnitude of the acceleration of the universe. The second problem is 
the so-called "coincidence" problem. At high redshifts, say just after 
inflation (z ~ 10 28 ), the radiation density was g ra d ~ 10 61 GeV 4 , while 
g\ <~ 10~ GeV 4 . This means that in order to have the correct amount 
of dark energy today, we need to fine-tune the initial conditions such 
that the ratio above was 10 110 with a very high precision. This also 
seems very unnatural. These considerations leaded to consider that 
the dark energy is time-dependent. As is known from inflation models, 
a simple way of obtaining these features is to consider a minimally 
coupled scalar field ip (termed quintessence). The latter is not neces- 
sarily the inflaton field, because we are at a much smaller energy scale, 
although the dynamical equations can be similar. Since the scalar field 
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evolves with time, the equation of state itself varies with redshift, and 
it constitutes a very important observational challenge to measure this 
dependence [6]. 

It is beyond the scope of our work to analyze the whole features of 
the quintessential models, but instead we elaborate on the suggestion 
of Bean and Magueijo [7], who suggested an important growth of seed 
PBHs when a quintessence scalar field ip is absorbed, dominating the 
accretion. The key new ingredient is the role played by the kinetic term 

. 2 

(^-) in the flux onto the PBHs, which is absent in the case of accretion 
of pure radiation. They further showed that the flux at the horizon 
is given by F v {t) = %f(t) (in natural units), and the potential term 
V(ip(t)) drives the behavior of the kinetic term through the dynamical 
equations describing conservation and the cosmic expansion. 

Type la supernovae data suggests that at z < 1 the dark energy is 
dominant [8]. It is sure that for z ^ 1, the dark energy was not domi- 
nant, because eventually matter becomes important at these redshifts 
and earlier, for z <~ 10 4 , radiation becomes more important. Then, we 
can conclude that, although we can not yet measure the dark energy 
contribution at z > 1, that contribution is subdominant but may be 
not negligible for the black hole growth problem from epochs earlier 
than present. We can view this in other form. The complete dynam- 
ics of quintessence plus radiation (and some black holes mixed as a 
secondary dark contribution) shows attractor solutions. The simplest 
models with quintessence show that it is very hard to have at, say, 
z ~ 10 and flDE « 1, while today Qde(z <~ 1) <~ 0.7. There are not 
"natural" solutions containing these features. 

We shall describe a system considering that the energy and pressure 
of the quintessential field dominate the black hole contribution, and one 
simplified model arises from there. 

A general model containing black holes + quintessence of the type 
studied here is described by 

S = J d 4 xy/^g(nR + L v + L matter ) ; (2) 
where the quintessence scalar field has pressure and energy-density 

•2 -2 

given by P(t) = \ — ^(v 5 ) an d g(t) = + V(ip) respectively. 

We assume that our spacetime is asymptotically flat, consisting 
of this quintessential field and black holes (without charge or angular 
momentum) absorbing the energy of this field and neglect all other 
contributions, here k = l/lQitG is the gravitational constant. The 
total does not need to be diagonal, even if the metric is diagonal 
in the beginning. The non-diagonal terms describe energy-flux at black 
hole horizons. 
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The dynamics of the system is given by 



dM 



CM 2 



■2 



growth 

(p + 3Hip + V' = ; 



(3) 
(4) 



H2 =(S)=^V + Spbh )-^-, (5) 

where a(t) denotes the scale factor of expansion/contraction of the 
universe and K = — 1,0, +1 describes the cosmic curvature. We 
shall evaluate physical effects in the flat model K = for simplic- 
ity, and assume that all the black holes started with the same mass, 
i.e. a single-mass scale initial function. We use natural units in 
which C — 277r/M p / 4 . Note that eq.(4) is the conservation law for 
the quintessential field in the form g + 3H(g + P) = 0. 

Since we have assumed from the scratch that PBHs do not domi- 
nate the dynamics of the universe g p bh « V(ip), analytical solutions 
for ip(t) and its derivatives can be obtained without additional com- 
plications, ft is obvious that, when this requirement breaks down, 
ip(t) becomes a very complex function of t and M(t) and a numerical 
analysis is necessary. 

Bean and Magueijo [7] have recently analyzed only one particular 
model of quintessence in which V(ip) — Xexp[\ip], implying ip 2 oc t~ 2 
within the slow-roll approximation. Within this model, in which the 
quintessence flux around the SMBH decreases with time, they claim 
that the black hole mass grows with growing time. However, the de- 
creasing quintessence flux around these objects (oc t~ 2 )implies that the 
mass gain term decreases accordingly (M cx M 2 t~ 2 ). In fact, we have 
proved that in the radiation-dominated era the growth of black holes 
is quenched when the background flux decreases with time as t~ 2 or 
faster (this argument is generalized in the present work for quintessence 
accretion, see also Ref.[9]). It is important to gain some feeling about 
the behavior of simple solutions first, and we start with one of the sim- 
plest possibilities, namely that the flux of quintessence onto the black 
holes remains constant. 



5 



2.1 Constant quintessence-flux through the hori- 
zon event 

Since we are interested in evaluating the effect of the quintessential 
field on the black hole growth, we may start solving the particular case 
in which the kinetic term (and hence the flux) F(t) oc ip 2 is constant. 
If we set ip = constant = B, cq.(4) yields 



V487ra(B) _ dV 



M, 



P i 



(7) 



with a(B) 



V487T-B 
M„i 



and a slow-roll dynamics of the quintessence 



field holds. Solving for V yields 



with a(B) 



V(ip) = V Q y/l + *(B)(tp-tp ) ; 

I 12tlB \1/2 



(8) 



The scalar field itself evolves as (p(t) = ip + v2B{t — to). Thus, 
the explicit solution for the evolution of the mass of the black hole is 



M(t) = 



Mi 



1 



M 



P l 



(t-U) 



(9) 



For t^S>ti, eq.(9) becomes 
M(t) = Mi x 



I - 2.<S / 10 '■■!,„ : ( — ) MT 



B 



(10) 



where ^ = (M t /M & ), r = (t/t ) and g c (t ) ~ lO^^gcm" 3 . 
It is clear that the asymptotic value of the mass is much larger than 
the initial mass M(t) ^> Mi if the denominator above satisfies 



< 



1 - 1 - 2.8 x 10 



< S 



(11) 



A/' 



with i5< 1. The PBH mass actually diverges at t a (Mi) — 54wBM . , 
numerically given by 



t a (Mi,B) ~ 3.4 x 10 



21 



"£»c(*o)~ 






B 


X 


M _ 



to ; 



(12) 



with to ~ 10 17 /io and g c (to) ~ 10 _29 gcm -3 . This of course does 
noi mean that a black hole will achieve an arbitrary large mass, but 
rather that before this happens, our approximations will break down 
and the condition ip = constant will not be satisfied. An analogous 
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exercise in which the quintessence flux decreases as F(t) oc t -1 still 
produces a (slow but reasonable) growth over the cosmic time. How- 
ever, for F(t) oc t~ n , with n > 1 the growth is quenched because the 
environment does not provide enough fuel as time goes by (see Section 

3). _ 

The next natural question is: how do the PBHs avoid the divergence 
of their masses in a "natural" way? For a generic quintessential field, 
but keeping the general slow-roll approximation for the dynamics of ip, 
we use the eqs.(3), (4) and (5) to obtain the mass of the black hole as 
a solution of a more complicated equation 

UM\ _ CM pl 2 M 2 (V'(p(t))+pf 

\dt) 48tt V(<p(t)) + g pbh (t) • 1 ' 

The non-trivial effects of the presence of other black holes are now 
explicitly seen. In the limit g p bh > V{<p) it is expected that the flux at 
the horizon decreases and the relation F(t) cx (p 2 ceases to be valid. On 
the other hand, the limit V(if) 3> g p bh leads to very simple solutions 
since this equation is separable (note that g p bh{t) = Mn p bh) and con- 
sists of growing solutions, holding for a finite time while V(ip) > g pb h- 



2.2 Evaluating the fractional BH-mass variation for 
quintessence flux. 

Let us evaluate the fractional change in the black hole mass (born with 
an initial mass Mi) when the quintessence field alone fuels the growth. 
Writing f2 v = — and noting that the contribution of the kinetic 

term is then quite generally we have £l v — 2 ^ . 

Therefore, g c Q v ~ \, whenever ^ip > V((p). When the kinetic 
term dominates, we have ^ = jtt-M 2 g c (t)tt v (t), where t = t(z). 

This equation is easy to solve. Integrating it from t(zi) until t(z = 0), 
and from Mi until M(t) yields 

Since the variation of the £l v is expected to be relatively small at low 
rcdshifts, and the quintessence term is known to amount ^0.7 today, 
we can always estimate the kinetic contribution of the quintessence as 
a function of the matter and radiation contributions. Neglecting as 
before the curvature term we may rewrite the eq.(14) above as 
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/AAA 2.77TM, f° , dt' 

{-w)~^yjj zgc{z) d? ; (15) 

which, after substituting the relations between z and t, in terms of 
H , becomes 



/AAA 



{zi) ~ MiH oZi ; (16) 



which is linear in Zi, as expected. The formula above is a lower 
bound estimate, because we know that Q v actually grows in that pe- 
riod. Then, it constitutes a minimal bound, when that variation of the 
latter is small. 

If we further parametrize the flux as F v = <fo(t/to) for a generic 
index n, we can obtain from the condition Q < 1 



¥>o < VQc(to) ~ W- 23 GeV 2 . (17) 

Supposing that the kinetic term dominates over the potential one, 
and inserting this upper bound back into eq.(16) we obtain 



dr J 



I < 1(T 2 VA ; (18) 



where n = and r = t/t - 

In these particular parametrized power-law models, n can be posi- 
tive, negative or null, and the results will differ for each case. Consider 
n > first. We can rise the following question: when will the net mass 
variation be relevant compared to the initial mass?. Let us suppose 
that we want to obtain a fractional growth of ~ 10 8 . Inverting the 
above relation, the cosmological time when this huge final mass will be 
achieved is 

t{}M) - io 3 °/( n+1 Vr 1/( " +1) to ; (i9) 

in other words, if the conditions fl v < 1 hold for future times and 
n is constant and positive, only in the very remote future the scalar 
field accretion will be able to accrete 10 8 times its initial mass. This 
future epoch of the universe will be coasting (provided tttotai = 1), and 
dominated by big black holes that feed from the quintessential energy 
in a very slow way. For n < 0, these conditions become much more 
severe, and we do not get any substantial growth unless very strong fine 
tuning had happened in the remote past. These conclusions are not 
likely to change for realistic (or more complete models) if we find the 
general background flux solutions given by the Friedmanns equations, 
since the fractional mass gain is 
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) 



10 



22 




(20) 



and /(t) is a consistent solution of the complete set of equations. 
This is because it is very difficult to cancel the tiny coefficient 10~ 22 , 
which arises quite generally in the consideration of the cosmological 
accretion. That is why we state that extreme BH growth can arise 
only in very peculiar situations, and why we believe that this is very 
difficult to achieve. 

One can also invert the reasoning above and evaluate the density 
of radiation/quintessence in remote past which would have produced 
(tT") 1 an d arrive to similar conclusions, namely that even given 
the most favorable conditions the quintessence would not be able to 
cause a huge growth of PBHs. 

3 Extreme fine-tuning for the growth of 
black holes immersed in a quintessence back- 
ground flux 

Using the critical mass (see [9] for discussion and references) as a tool 
let us show now that any flux of the form F v oc t~ n , with n > 1 does 
not lead to a substantial growth for any black hole, unless we require a 
strong fine-tuning between the formation-time and the initial mass of 
the black hole (with all other parameters kept fixed). This functional 
parametrization is specially suitable for most applications and permits 
to make a quick contact with standard cosmological evolutions. 
The dynamical equation in the accretion regime is 




(21) 



where C, and i* are constants. Plain integration yields 




(22) 



with r) = (t/t*) and E 1 * = t* 
The l.h.s. is just 



-l 



1 



1 



(23) 



Mi 



M{rj) 



and the r.h.s. is 
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CM pl 4 
(n - l)E. 

Then, the solution M{r}) is 



n n-l 



M{n) 



1 

Mi 



1 - 



CM z M pl 4 



,71-1 



n-1 



Considering n > 1 we obtain for asymptotic times, i.e r\ 3> r\i 
Mi,) = 

(n-l)E, 



(24) 



(25) 



(26) 



It is clear that the denominator is much smaller than 1 (and pos- 
itive) at the expense of a very strong fine tuning between Mi and t%. 
Let us evaluate a definite numerical example. Suppose that the growth 
factor must be of order 10 s as before, i.e the denominator 



1 - 



CM pl 4 Mi 
(n - 71"- 1 



= 10" 



(27) 



Then, the formation time must have been very fine-tuned into 
the function above (in one part by several millions) in order to get 
a huge growth factor. If this did not happen then either the growth 
factor falls off quickly to 1 or the solution is unphysical. Note that this 
strong conclusion is very general, and it is independent of the origin 
of the flux through the black hole horizon, i.e. the result is valid for 
the radiation-dominated era, in which the flux was cg rac i(t) cx t~ 2 or a 
quintessence-dominated era when tp 2 cx t~ 2 (see Ref.[9]). 



4 Further constraints on PBH growth from 
the Generalized Second Law and the Holo- 
graphic Bound 

The irreversible process of quintessence disappearance into a black hole 
is just a particular example of the Generalized Second Law of thermo- 
dynamics (hereafter GSL) formulated in the '70s by Bekenstein [10] 
and other following works [11, 12]. 

The GSL is the statement that the total entropy (black holes plus 
any other fields) never decreases with time 

AStotal = AS bh + AS matt er + AS^ + AS rad i at i on > ; (28) 
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where 



S bh (N,M)~W 77 N(M/M Q f ; (29) 

counts for the entropy of ./V-black holes, S ma tter and S ra diation is 
the normal matter and radiation entropies, and S v is the entropy of the 
scalar field. A step forward towards the use of entropy as a tool in cos- 
mology is the so-called Holographic Principle (HP). The Holographic 
Principle [11, 12, 13], is the statement that all entropy/information 
contained in the physical universe can be attributed to its boundary. 
In numerical terms, the upper bound to total entropy (known as the 
Holographic Bound) is 

s max (t) < ^jt ; (30) 

where A(t) is the boundary area of the physical system at t. If wc 
assume that the HP can be applied to the cosmological models, then 
it is natural to choose the cosmological particle horizon to bound the 
area and write the maximum entropy as 

Sft P (t) ~8x 10 121 (tAo) 2 • (31) 

If the inequality Sbh(N,M) < Sh P (t) has to be satisfied, the accre- 
tion of quintessence should adjust in such a way that a black hole with 
initial mass Mj may grow to values Mf > 1O 6 M , without violating 
the Holographic Bound. Therefore, an upper bound to the local flux 
due to this global requirement is obtained. 

It is clear that this energy input strongly depends on the initial 
mass of the PBH Mi and needs to be very large if the black hole was 
initially small. In addition, if the accretion is very high, the black hole 
that was initially below the Holographic Bound will blow that bound 
at some point. Then, to keep these black holes below the Holographic 
Bound at any time, the constraint 

S S mbh{t) < Shp{t) ; (32) 

must be satisfied, together with S sm bh(ti) < Sh P {ti) at the forma- 
tion time U. When both conditions are required to hold, an absolute 
upper bound to the flux onto the N accreting black holes can be es- 

, or numerically 



M 

tablished. It comes out to be F v < -# 



10 aa t 
CNM 3 (t) 



(t/s)h 2 
{N/lQ ll ){M/M^' : 



F v {t) < 4.8 x l0- 5 7T77T ^Z^__ 5cm -2 s -i . (33) 
with M* - 10 8 M m . 
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If the total number N of black holes increases, the flux F v (t) would 
be correspondingly smaller, in order to keep all black holes growing at 
moderate rate and the total entropy smaller than Ahp/A. 

A general and rigorous proof that the Holographic Bound can be ap- 
plied to the cosmological models is lacking now, but if we proceed with 
these ideas we can get some insight about the local and global proper- 
ties of black holes in realistic models of the universe, although the full 
meaning of eq.(44) is obscure yet. Some problems of these methods in 
cosmology are discussed in Refs.[12] and [13]. Upper bounds to black 
hole abundances have been obtained using these ideas [14]. 

Quite regardless of the use of the Holographic Bound, we expect 
that when the mass-energy of black holes is comparable to term V(ip) 
in the Fricdmann equation, the system response would be to deplete 
the quintessence flux onto the black hole horizon by decreasing Lp 2 . A 
complete study of these features will be published elsewhere [15]. 



5 Conclusions 



We have shown in this work that, generally speaking, the quintessence 
flux ip 2 must decrease slower than t~ 2 for PBHs to grow at all, and 
to stay constant or increase in time for substantial accretion to occur. 
This may be relevant for some PBHs achieving the SMBH condition 
in a short time, but only at the expense of fine-tuning in the parame- 
ters. Several options for the quintessence potential V(ip(t)) have been 
proposed in the literature, some of which can lead to high flux to feed 
PBHs (i.e. leading to large positive M) for long times. One analytical 
solution of this class leading to ip 2 = constant has been worked out as 
an example. 

If the flux onto the black holes is constant, all them initially above 
the critical mass will grow until the term qpbh approaches the energy- 
density in the cosmic background. From this point, these black holes 
would not keep growing at high rates, and the more complicated cal- 
culations have not been performed for this regime. 

In Section 3 we showed that if the flux behaves as F{t) oc t~ 2 , 
we can prove that all solutions with asymptotic growth (M(t) ^> Mj) 
depend of an initial fine-tuning between the initial masses M, and their 
formation at t = U. Other fluxes that decrease more slowly (or even 
increase) than t~ 2 may lead to very big black holes but only at t ^> U, 
a condition that is not favorable for quintessence accretion which is 
unimportant at high redshifts and thus does not have enough time to 
happen. 
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In Section 4 we have invoked Holographic Bound arguments to show 
that a bound on the quintessence flux can be obtained quite generally, 
although these arguments would only hold if a rigorous proof of the 
existence of a bound is obtained. In spite of the theoretical efforts, this 
has not been achieved and therefore the limits must be considered as 
speculative. 

It should be kept in mind that models others than the one discussed 
here can be constructed to produce a population of SMBHs starting 
from seed PBHs. For example, accretion in a brane-world high-energy 
phase has been recently studied by Guedens, Clancy and Liddle [16] 
and Majumdar [17], which show that a substantial growth in which 
M cx M/t is allowed. It may be possible to arrive to the end of 
the high-energy phase with very massive black holes, although the full 
consequences of this type of scenarios are yet to be explored. 
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7 Appendix: the critical mass in radia- 
tion+quintessence environments 

A black hole surrounded by radiation/quintessence would be either 
absorbing energy or evaporating, according to the sign of the differ- 
ence Tpbh — T environment , just like any other thermodynamical system. 
Consideration of the complete differential equation describing the black 
hole mass evolution 



is needed (where F to tai (t) represents the total flux of energy through 
the black hole horizon (that is, the addition of fluxes for any form of 
energy that falls through the black hole) at cosmic time t). For a 
given Ftotai the condition ^M. = defines a parameter (the critical 
mass M c (t)) separating the absorbing and evaporating regimes. The 
most interesting cases arise when the flux is dominated by radiation or 
quintessence. We shall discuss the cases separately to highlight some 
interesting differences. 

7.1 Radiation-dominated environment 

The case of pure radiation Ftotai = F ra d(t) = cg ra d(t) has been ana- 
lyzed in [4]. By solving the algebraic equation ^ = 0, the critical 
mass is found to be 




(34) 



M c {T rad ) 



6 1 26 ff 



(35) 



(t) (T rad /T ) 
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where 6 = M pl (A/27ira*) 1/4 , g rad (T) = a*T 4 and T ~ 3K is the 
present temperature of the CMBR. 

Then, if we ignore other possible sources of energy-matter to fuel 
their growth, we conclude that any black hole with mass > 10 26 g is 
today in the absorption regime. The critical mass has evolved with 
time, it was M c ~ 10 16 (/ when the universe was t ~ Is old. Therefore, 
any PBH with initial mass greater than the latter value formed within 
the radiation-dominated era, was initially above the equilibrium con- 
dition, and the Hawking evaporation was irrelevant for it. Apart from 
fine details of the absorption process, the cross section is expected to 
be valid quite generally (see Refs. [4] and [9] for a full discussion). 



7.2 Quintessence-dominated environment 



If we consider the quintessential field as the dominating fuel, the bal- 
ance equation for the black hole mass ^M- = becomes 



(dM 
\~dT 



total 



M 2 



+ ■ 



CM 2 



(36) 



In order to simplify the expressions, let us define the dimensionless 
ratio A = oM^ <~ 1.1 xlO -3 . Thus, the critical mass for this case is 
expressed as 



0.06M p , 

*.<«) = ^rjr 



5.6 x 10 24 9 



2g c (t ) 



-1/4 



(37) 



where we have defined the dimensionless function E(t) = 2 ^ 4. ■ 
The expression shows that if the quintessential density is of order 
of the critical density today, the corresponding critical mass is about 
10 25 g. This is not very different than the value corresponding to pure 
radiation dominance, a fact directly related to the assumption V ~ 

Let us show how the constant flux case described in Section 2 looks 
like in terms of the critical mass. The critical mass is a constant given 

by 



M c (t) - 5.6 x 10 24 g 



B 



-1/4 



(38) 



Therefore, only initial masses greater than this value can be de- 
scribed by growing models with constant quintessence through the 
horizon. Fig. 2 displays the behavior of black holes in this situation 
according to their initial mass, see caption for details. More general 
models can be studied. For instance, let us briefly discuss the case of 
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a power-law flux F(t) = M pl {t/t^) n with t* a reference time. In the 
case n < 0, we obtain for the critical mass 

M c (t) ~ (<7/277r) 1/4 M p *(W H/4 ; (39) 

and this critical mass increases with increasing time. Since the 
ratio of the absorption term to the evaporation term is just 

iU? = (M/M C ) 4 ; (40) 



(valid for any form of the accreting flux) then, if we follow the 
evolution of a PBH mass whose initial mass Mj was slightly larger than 
M c (t — ti), then we can ignore the Hawking term on its evolution until 
the PBH mass crosses the critical mass curve. The error introduced 
to obtain the crossing time t c is small. Thus we can describe the PBH 
mass evolution as 



1 1 27ttU 



Mi M c (t c ) (n - 1) 



1 1 



n 1 _n— 1 



(41) 



where rji = and r\ c = (t c /t*). Since t c ^> U in general, the 

second term is always negligible. 

EDITOR: PLEASE PUT FIG.l HERE!!!!!!!! 

Solving eq.(22)for the crossing time yields 



4/rt / \ -l/n 



l-^(^) 1 - 

n — 1 



-4/n 



(42) 

Then, in a finite time, this PBH will achieve an instantaneous equi- 
librium at the critical mass crossing, and immediately proceed into the 
evaporating region. 

In the case n > 0, the critical mass is 

M c {t) ~ M pl (a/277r) 1/4 (t/uy M/4 . (43) 

And this function decreases with increasing time. Suppose that a 
PBH has been formed with initial mass smaller than the critical mass 
at ti. Then, for a long time this PBH will loss energy until its complete 
evaporation (if we ignore other quantum corrections likely to arise near 
the Planck scale). However, if its mass reached the critical mass at 
value M c (t c ) > M p i, then the black hole achieved an instantaneous 
equilibrium at t c (see below) and it will proceed into the absorbing 
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region, from which it will be growing increasingly faster (Fig. 3). We 
note that at M = M c we have M(M = M c ,t c ) cx M 2 (t/t*)" _1 > 0, 
and it will be increasing for n > 1. Now, if a PBH had already formed 
below the critical mass at ti, and is not much smaller than it we can 
solve the mass differential equation to find the crossing time 



This solution is numerically valid only when i c (Mj) <C t evap ~ 
t (Mi/M haw f. 

From t c (Mi), this black hole will be accreting and its mass will 
be growing. These considerations illustrate the role of the critical 
mass for the evaluation of the fate of PBHs. Any statement about 
growth/evaporation is necessarily time-dependent and must be care- 
fully done to preclude a misunderstanding of the physical situation. 

EDITOR: PLEASE PUT FIG.2 HERE!!!!!!!! 

In a completely generic case F to tai = F(t) we obtain 
2A 



M{M, t) = 



M^ 2CBM 



M{M,t)+CM z {t)F{t) ■ (45) 



and the fate of these black holes will depend on the sign of F(t) 
as expected. In general, M(M — M c , t) is not zero, and a similar 
conclusion holds for higher derivatives as well. Therefore, we conclude 
that very special conditions have to be fulfilled for black holes to evolve 
in equilibrium with the environment, that is, without losing or gaining 
mass. 

7.3 Quintessence/radiation environment 

If both radiation and quintessence are considered, then the critical 
mass becomes 



M c {t) 



10 



1 + — 



Brad 



-1/4 



(46) 



(Trad/To) 

to be solved together with eqs.(4-6). In general, the behavior of 
Brad would not be cx t~ 2 , which happens only when the radiation 
dominates all forms of energy. In this case, the first derivative of 
M(M, t) at M = M c is given by 

M(M = M c , t) = CM 2 (ATT 3 + (pip) ; (47) 



17 



and this term is in general not zero unless pip = —4TT 3 . Since 
radiation and quintessence would have to conspire for that, equilibrium 
of PBHs and the environment would not be possible. In scenarios 
where V(ip) >> g ra d PBHs can stay in equilibrium, actually one of 
these is the constant flux model already discussed p 2 = 2B = cte. 

Now, we can evaluate the redshift dependence in these terms. In 
any way, g ra d(z) oc (1 + z) 4 , regardless of whether the quintessence 
dominates. We can cast the eq.(47) above as 

10 26 o 

M c (z) 9 - — ^ • (48) 



( l + z ) + eJ,(0) (df 



It is easy to conclude that if -jj^yc- — j^J ^Zz ^ z 2 oc (1 + z) ra 4 

the following behaviors are possible. If m > 4 the quintessence will 
determine the critical mass even at high redshifts. For m < 4, the 
opposite case happens, and there will be a redshift z; from which the 
radiation will determine the critical mass. It will be true at (1 + zi) ~ 

"I 1/4 

, evaluated at zi also. It depends (as it should be) on the 



it)* 2 

time derivative of the quintessential field, and the complete solution 
must be solved when we choose some form of V(ip). It is important 
note that when the term V(<p) is relevant, the accretion into black holes 
can be null in some cases even when Q, v is big. We can write Q, v = 

' f \ Then, when the first term is important, the quintessence 



Be 

dominates and the black hole accretion (from it) is important also, but 
the inverse is not true. We can choose other models with V(tp) leading 
to \(p 2 = 0, and the quintessence dominates without some secondary 
effect onto black holes with any mass. The third point now is. When 
the form of V(ip) leads to a non-null kinetic term, the magnitude of 
this term depends also on the magnitude of V(<p) (obtained by solving 
the dynamics) and the accretion onto black holes will depends on the 
parameters that enter into V(tp)- definition. Finally, the fourth point 

is that the z-dependence in Vl v (z) = 2 ^ ^^j^^ reflects in each 
term. Today we know that fl v (z <~ 0) <~ 0.7, but the attractor features 
of the solutions indicate that when z > 1, the magnitude of the terms 
above were not so much below the unity also, although the exactly 
dependence and values are not known yet. In terms of time, from 
z = 10 to z = there was enough time to accretion into black holes, if 
the others conditions were satisfied together also. Moreover, there are 
indications that the quintessence varies with z in a very slow way if at 
all, see [18]. 
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7.4 The behavior of M(M = M c , t) 



An interesting related question is whether PBHs formed with Mj = M c 
can stay in equilibrium without evaporating or gain. We give a brief 
assessment of the general features of the solutions M(t). Given the 
general form of M(M,t) we may express M(M c ,t) in the following 
way 

M(M c ,t) = CM 2 (t)F(t) ; (49) 

in terms of the mass and the behavior of the flux at its horizon. 
Since M = M c , then it will be simplified to CM 2 (t)F(t) and it is 
enough know the evolution of the critical mass and flux. The true 
behavior of each black hole mass is very complicated even in the simply 
cases above, but through the critical mass and its derivatives, we may 
infer the black hole mass at asymptotic times. Then, the algebraic 
sign and derivatives of M(M c ,t) will determine the behavior of the 
solutions. 

Applying these definitions to the constant flux case we obtain M(M - 
M c , t) = on the critical mass. The first derivative of M(M, t) is just 

ft { ~ -m + CM 2 B), rewritten as 



M{M,t) = 



OA 

+ 2CBM 

M 3 



M(M, t) . (50) 



Since M(M = M c , t) = 0, equilibrium is possible on the critical 
mass and the black holes will retain their initial masses. This result is 
valid as long as we ignore the radiation compared to the quintessence 
field. Note that all the higher derivatives of M(M,t) are null at M c 
only if the flux through the horizon is constant. 

For the example mentioned in Section 3 of a power-law flux F(t) — 
M pl(t/t*) n - The critical mass is M c (t) = pl f i{ with a = (A/C) 1/4: 
and the second derivative M reads 

M(M c ,t) = —^(tjuf l2 - Y . (51) 

From eq.(31) above, it is clear that we have a branching of the 
behavior for the solutions. The first branch corresponds to n > 2 and 
in the limit t — > oo we obtain M(M c ,t — > oo) = oo. On the other 
hand, n < 2 yields M(M C , t — > oo) = These cases are physically very 
different, because the function M(M, t) is related to the energy that 
crosses the horizon event, and its derivative is sensitive to changes in 
the asymptotic limit. The cases < n < 2 correspond to solutions 
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showing moderate growth and the case n > 2 to very fast growth of 
the black holes. The cases n < are the ones for which the critical 
mass grows, the behavior of black holes in this latter case is determined 
by the initial condition Mj. 
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Fig.l. The critical mass for a constant quintessence flux. The critical mass is 
represented by the horizontal line labelled in the graph. A black hole may either 
grow or evaporate according to the value of its initial mass Mi above or below the 
critical mass respectively. Instantaneous equilibrium of a black hole not growing 
nor evaporating is possible if its mass is precisely the critical value, as discussed in 
the text. Note that growing is expected to saturate when the primordial black hole 
density approaches the quintessence field density. 

Fig. 2. The same as in Fig. 1 for a power-law quintessence flux of the form 
F(t) = Mp^t/t*) 71 , with n > 0. The black holes suffer a delay of their evapora- 
tion and may start to grow if they manage to cross the critical mass curve which 
decreases with increasing time. This case illustrates how different acceptable mod- 
els of quintessence have to be analyzed individually to address the fate of PBHs 
possibly formed at early times . 

Figure captions. 
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